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Abstract 

The existence of a nuclear electric octupole moment (EOM) requires both 
parity and time invariance violation. The EOMs of odd Z nuclei that are 
induced by a particular T- and P-odd interaction are calculated. We compare 
such octupole moments with the collective EOMs that can occur in nuclei 
having a static octupole deformation. A nuclear EOM can induce a parity and 
time invariance violating atomic electric dipole moment, and the magnitude of 
this effect is calculated. The contribution of a nuclear EOM to such a dipole 
moment is found, in most cases, to be smaller than that of other mechanisms 
of atomic electric dipole moment production. 
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I. INTRODUCTION 



At present, time invariance violation has only been observed indirectly, in the CP- 
violating decay of neutral K-mesons M. Parity and time invariance violating nuclear 
and atomic multipole moments, such as the magnetic monopole, electric dipole, magnetic 
quadrupole, and electric octupole moments, are interesting because, if discovered, they would 
provide further proof of time invariance violation. Even the limits on these moments provide 
important tests of different models of CP-violation. Parity and time invariance violating nu- 
clear moments induced by T- and P-odd nuclear forces have been discussed, e.g., in Refs. 

M- 

In this paper we consider the electric octupole moments of nuclei. We calculate the 
EOMs of odd Z nuclei that are induced by a T- and P-odd interaction between the unpaired 
proton and the nuclear core. Note that nuclei with unpaired neutrons can have an EOM of 
comparable magnitude due to the polarization of the nuclear core by the T-, P-odd field of 
the external neutron. A similar mechanism for the electric dipole and Schiff moments was 
considered in Ref. ||. 

We give values of the collective electric octupole moments of nuclei having a static oc- 
tupole deformation, using the calculations done in Refs. |||| . We also present a calculation 
of the atomic electric dipole moment (EDM) that would be induced by a nuclear EOM. 
Finally, we discuss a possible enhancement mechanism for magnetic quadrupole moments 
(MQMs) in nuclei with octupole deformation. 

In the appendix we present simple estimates of the relative sizes of the contributions of 
various nuclear moments to the atomic EDM. The magnitude of the nuclear EOM is com- 
parable to that of the Schiff moment. However, the contribution of the EOM to the atomic 
EDM is smaller than the contribution of the Schiff moment (and the magnetic quadrupole 
moment) since the EOM interacts with higher angular momentum electron states, whose 
wave functions are suppressed near the nucleus. 
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II. THE HAMILTONIAN OF THE T-, P-ODD NUCLEON-NUCLEUS 
INTERACTION AND THE RESULTING NUCLEON WAVE FUNCTION 



For a heavy nucleus, the T- and P-odd interaction between a nonrelativistic unpaired 
nucleon and the nuclear core can be described by the following effective Hamiltonian (see, 
e.g., Refs. 

where er is twice the spin operator for this nucleon, p is the density of the nuclear core, G = 
1.0 x 10~ 5 /m p 2 is the Fermi constant, m is the mass of the nucleon and r] is a dimensionless 
constant that describes the strength of the interaction. In this paper we deal with odd Z 
nuclei, and so the nucleon involved is the unpaired proton. 

Let U be the strong nuclear potential of the core that the unpaired proton moves in. 
The range of the strong nucleon-nucleon interaction is small. This means that the potential 
U(r) and the nuclear density p(r) will be similar in shape. In fact, we assume that they are 
approximately proportional: U(t)/U(Q) ~ p(r)/p(0). This allows us to rewrite Eq. ([!]) in a 
form that makes it easy to find the perturbed wave function due to Htp (we use the method 
of Ref. §). We obtain 

H TP « £o- ■ VU, (2) 

where 

< - "^LM = 2 x 10 ~ 21 " cm - (3) 

The total potential that the unpaired proton experiences is 

U = U + H T pmU(T) + Ztr-VURiU(T + £(r). (4) 

As a result, if ^>(r) is the proton's wave function when only U(r) is present, the perturbed 
wave function will be 

^(r) ^^(r + ^cr) ^ ip(r) + £(T ■ Vip(r). (5) 



3 



III. THE ELECTRIC OCTUPOLE MOMENT OF A NUCLEUS WITH AN 

UNPAIRED PROTON 



In this section we calculate the electric octupole moment (EOM) of a nucleus with an 
unpaired proton using the perturbed wave function obtained above. The electric octupole 
moment can be written as [see Eq. (AH))] 



O ijk = (ij>\dijk\i(>) = e J ^ (r)[rirjr k - \r 2 (ri5 jk + Tjb ik + r fc <%)]^(r) d 3 r. (6) 

Note that we use e > 0. Substituting -0(r) = ^>(r) + C,cr m ^- (from Eq. (|)) into the above 
integral and expanding gives 

O ijk = 2e£ J -^ r o m \r i r j r k - \r 2 {ri5 jk + rjS ik + r^)}^ d 3 r. (7) 

Notice that we have discarded the £ 2 term and that the term not containing £ vanishes, as 
it is an integral of an odd function of r. Applying integration by parts to the above integral 
gives 

O ijk = -\e^\h{fifja k + hhoj + rjf k &i) - r 2 (&i5 jk + ofi ik + cr k 5ij) 

-2r m a m (fi5 jk + fj8 ih + r fc <%)|^)- (8) 

We can also write the octupole moment tensor in another form. Since it is a symmetric, 
irreducible (traceless) third rank tensor and the nuclear angular momentum I is the only 
quantity which defines a direction in the system, the EOM tensor must be in the form 
of the most general symmetric, irreducible third rank tensor that can be formed from the 
components of 1. That is 

O ijk = (^|4i*#>, (9) 



where 



Oij k — A[liljlk + Ijlkh + Ikhlj + Ikljh + Ijhlk + hlklj 

- 6I(I+ 5 1) - 2 (i i 5 jh + i 3 5 lk + / fe( y] (io) 
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and A is some constant. The factor of —[61(1 + 1) — 2]/5 follows from the requirement of 
tracelessness (Ouj = = Oju = 0). The quantity which is usually referred to as the 
octupole moment is O, which is the O zzz component for the nuclear state having angular 
momentum projection I z = I. This is the quantity which we will calculate. We can write A 
in terms of O using Eqs. (^) and fllPP, with i = j = k = z. This gives 

Oijk 6/(7— 1)(27— 1) \ i j k ^j^k^i ^k^i^j ^k^j^i ^j^i^k ^i^k^j 

_ 6 /(/+i)- 2 ( ^ fc + + (n) 

We now have two expressions for Oij k : Eq. (||) and Eq. ([11]) (via Eq. (Q)). We calculate O by 
operating on both of these equations with on the left and I k on the right and equating 
the results. 

To evaluate the results of these operations, the following commutation relations are re- 
quired: [fi,crj] = 0, [ii,f-j] = isijkh, and = ie ijk a k . (Note that I { = k + &i/2.) 
Other useful relations are [I^Ijfj] = [I^Ijaj] = [I^Ujfj] = and e^AiAj = \sij k AiAj — 
^EijkAjAi = ^Eijk[Ai, Aj], for an operator A^. The commutation relations are used to rear- 
range the operators so that pairs with the same indices are adjacent. 

Applying the operators to Eq. (H) gives 

$\iiijd ijh i k \ij)) = -ge£(VI ^iJjififjCTk + fifkCTj + rjr k °i)h ~ ljjr 2 (ai5 jk + &jb ik + <x fc c%)4 

iT rn (J rn \T 

iSjk + rjSik + h5ij)h\4>) 
= -K(r 2 ){(/, /„ • n)[|(i • cr) - J 2 ](«x ■ n)|J, I z , I) 

+</, I t ,l\l(I-&)- 3/ 2 (I ■&) + !- 2/ 2 |J, I g , /)}, (12) 



where fi = f/r. Here the facts that I • n = (1 + <r/2) ■ n = (& ■ n)/2 and (<x ■ n) = 1 JTT 
were used. The radial and angular parts of the operators have been separated in Eq. (^2j); 
(r 2 ) is the expectation value of r 2 and \I, I z , I) is the angular part of The second matrix 
element in Eq. (|T2J) is 

(/, I z ,l\l(I-&)- 3/ 2 (i • <r) + 1 - 2I 2 \I, I g , 1} = [f - 3/(7 + l)](i - k) + 1 - 21(1 + 1), 

(13) 



where we have used the fact that 

(I • &)\I, I z , l) = [1(1 + 1) - /(/ + 1) + 3/4] | J, I z , I) 

= (1/2- k)\I,I z , I), (14) 

with 

k= (I +l/2)(-l) I+1/2 - 1 . (15) 
To evaluate the first matrix element in Eq. (|T2"D, we use the following identity: 

(&-h)\I,I z ,l) = -\I,lJ), (16) 
where I = 21 — I. We then have 

(/, h, l\(a ■ n) [f (i ■ &) - I 2 }(& ■ ft) | J, I„ 1} = (/, I z , T\\(l ■ &) - P\I, J„ I) 
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k)- 1(1+1). (17) 



(Note that (I ■ a)\I, I z , I) = (1/2 + k)\I, 7 2 , /).) Using these results in Eq. flTJ) gives 

<^A^W> = -fe£(r 2 )( K - |)(J - |)(J + |). (18) 

Now we find another form of the left hand side of the above equation by operating on Eq. 
(0). We have 

$\iiijdij k ik\i)) = 6 /(j_i)(2/-i) $\iJj(iJjh + ijhh + iJJj + iJjh + ijhh + iJJj)h 

- 6I ^- 2 i i i j (i i 8 jk + i 3 5 lk + y^iS). (19) 

This can be evaluated using the equation [ij, Ij] = iSij k i k to bring operators with the same 
indices together and Eij k lilj = \eij k Eij V l v = il k . The result is 

(tiliiijdijJS) = 0(1 + 1)(I + 3/2) (J + 2). (20) 

Equating this result with Eq. ([0|) gives the following result for the octupole moment: 

-3(«-3/2)(/-l/2) 2 
Usmg " 5(/ + l)(/ + 2) (r;e4 ' [Zi) 
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(The subscript refers to the fact that this octupole moment is due to a single particle, as 
opposed to a collective octupole moment.) The expectation value for r 2 can be approximated 
as § 

(r 2 > » |r 2 A 2 / 3 , (22) 

to an accuracy of about 10%, where r = 1.1 fm and A is the mass number of the nucleus. 
Using the value of r and Eqs. (H) and fll5f) gives 



for / = / + 1/2 



C sinE « 8.7 x 10- 9 A 2 /V (fm)3 x { 1+1 . (23) 

1 ;/-i/2)(/-i) 

(J+l)(J+2) 



(^^(/-D for I = Z- 1/2 



Observe that C sing = for / = 1/2. This is to be expected because is a third rank tensor, 
and so applying the triangle rule for the addition of angular momenta to Eq. (§) gives the 
result that nuclei having angular momentum less than 3/2 cannot have an octupole moment. 
Values of S ing for various nuclei are given in table 0, in terms of the parameter rj. 

IV. COLLECTIVE ELECTRIC OCTUPOLE MOMENTS IN NUCLEI WITH 
STATIC OCTUPOLE DEFORMATION 

In Refs. JSIU, a mechanism was suggested by which parity and time invariance violating 
interactions can produce collective T- and P-odd multipole moments in even-odd nuclei 
having a static octupole deformation [i.e. electric octupole moments in their intrinsic (or 
body-fixed) reference frames]. Such a deformation has been shown to exist for nuclei in 



the Ra-Th and Ba-Sm regions (for a review see, e.g., fl!2|). A similar mechanism, for 
the enhancement of the intrinsic electron EDM and other T-, P-odd interactions in polar 
molecules, was suggested in Ref . |13| . Below, we explain the mechanism by which a collective 



nuclear EOM can be produced and provide values of this EOM for various nuclei. 
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A. Parity and time invariance and octupole moments in the laboratory frame 



An electric octupole moment can exist in the nucleus's intrinsic frame without parity or 
time invariance violation. Yet if parity and time invariance hold, the expectation value of 
the octupole moment in the laboratory reference frame will be zero. 

Consider \IMK) and \IM — K), which are two almost degenerate states of the nucleus 
in the laboratory frame. / is the angular momentum of the nucleus, M is its projection onto 
the z-axis, and K is its projection onto the axis of symmetry of the deformed nuclear core 
(the z'-axis). (I is the sum of the unpaired nucleon's angular momentum, j and the nuclear 
core's orbital angular momentum, R.) These states can be written in terms of intrinsic 
states as 



\IM ±K) = y/^D^te, 6, 0)^(r')Xcore, (24) 
where £^±#(0,0,0) is a Wigner D-function (see, e.g., |I4|,[15| ), Xcore is the wave function 



of the nuclear core in the intrinsic frame, and iP±k{ t ') is the wave function of the unpaired 
nucleon in the intrinsic frame, with a z' angular momentum projection of ±K. (In the 
intrinsic frame, the nuclear axis plays the role of the usual "2- axis" in Quantum Mechanics.) 
Note that j and I have the same z' projection. 

\IMK) and \IM — K) do not have good parity, as K changes sign under a parity 
transformation. However, the following states do, and they form a parity doublet: 

il> ± = -^=(\IMK)±\IM-K)). (25) 
v2 

For these good parity states (V^II ■ n|^ ± ) = because K and —K have equal probabilities 
and this means that there is no average orientation of the nuclear axis in the laboratory 
frame ((^ ± |n|'?/' ± ) = 0). This is a consequence of time invariance and parity conservation 
since the correlation I • n is T-, P-odd. As a result of ('?/> ± |n|'?/> ± ) = 0, the mean value of the 
octupole moment (whose orientation is determined by the direction of the nuclear axis) is 
zero in the laboratory frame. 



S 



Now, a T- and P-odd interaction, Htp will mix the members of the parity doublet 
and The admixed wave function of the predominantly positive parity member of the 
doublet will be ip = ip + + atp~ or 

V> = -7=[(1 + a)\IMK) + (1 - a)\IM - K)], (26) 
v2 



where a is a mixing coefficient: 



a - (27) 



— £L is the energy splitting between the members of the parity doublet. The interaction 
Htp is given by Eq. (H). (A similar expression can be obtained for the predominantly 
negative parity member of the doublet.) This mixing yields, on average, an orientation of 
the nuclear axis along the direction of the angular momentum: 

(V>|I-n|V) = (ip\K\ip) = 2aK, (28) 
and this means that the octupole moment need no longer vanish in the laboratory frame. 

B. The magnitude of the collective octupole moment 

The collective EOM in the laboratory frame was derived in [§,0, with the following 
result: 

^ 4 /(/-!)(/ -1/2) 
a ° U ~ 5(/+l)(/ + 2)(/ + 3/2) a ° 3 ' intr - (29) 

Once again, observe that O co ii = for / = 1/2. 0^i ntr refers to the octupole moment in the 
intrinsic frame (O zzz = §03 ) and is given by |16|JI7| : 



3 , intr = eZR 3 -^=((3 3 + ~W -fcfc + -^/3 3 /3 4 + •••), (30) 

2V77T O V 7T llV 71 " 

where R = r A 1 ^ (r = 1.1 fm). /3 2 , fa, and /3 4 are parameters that describe the nuclear 
deformation; the surface of a deformed nucleus is 
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fl = fi [l+E^o(M)]. (31) 



1=1 



We will first present an order of magnitude estimate of a. K = I • n and Htp are 
both T-, P-odd pseudoscalars. Therefore, (iP+k\Htp\iP+k) oc K and so {iP-k\Htp\iP-k) = 
— (iP + k\Htp\?P+k) (this fact can be easily supported by model calculations). Using this fact 
and Eqs. and (|25|) we get (^~|iJTp|V ;+ ) = (V'+id-^rplV'+ii")- If V ; +ic were a good parity 
state this matrix element would be zero. However, due to the perturbation caused by the 
static octupole deformation of the nucleus (V3), it is a combination of the opposite parity 
spherical orbitals ipi t +K and ip2,+K (e.g., P3/2 and d 3 / 2 ): 

tp+K = + l4>2,+K, (32) 

7 = F F ' ^33) 

i/)i )+ k = Ri(r')Q jt i t+K (9' ',</>'), 

i>2, + K = R 2 (r')n jX+K (e', 0') = -R 2 {r')(cr ■ n')n^ +K (0', 0'), (34) 

where I = 2j — I. (Of course there will also be an admixture of other opposite parity states 
having different values of j. We neglect these states for simplicity.) Therefore, we have 

{ip +K \H TP \ip +K ) (ip li+K \H TP \ip 2 ,+K) 
a = E + -E_ = 27 E—El • (35) 

To estimate 7 we must first derive the form of V3. As in Sec. |J, let U be the strong 
nuclear potential that the unpaired nucleon moves in. V3 is the difference between the 
potentials with the octupole deformation present, ^7 pre s and absent, £7 a bs- Once again we 
have U(r') ps U(Q)p(r')/ p(0). We also make the approximation that p(r')/p(0) ~ 0{r' — R), 
where R is the nuclear radius and 6(x) = 1 for x < and 9(x) = for x > 0. For an 
octupole deformed nucleus we have R = Ro(l + PsY 30 ) and so 
V 3 = U pms - U ahs w U(0)[e(r' -Ro- (3 3 Y 30 R ) - 6{r' - R )} w U(0)R (3 3 5(r' - R )Y 30 , 

(36) 
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where we have expanded 6 in a Taylor series, using 4^ = —5(x). Using Eq. (^) we then 
have 
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U(0)f3 3 R 1 (R )R 2 (R Q )R 3 [ ifecA dn\E l - E^ 1 



Pb, (37) 



where we have used R l (R )R 2 (R ) ~ l-4/i? 3 0, 1^(0) | ~ 50 MeV, \E X - E 2 \ « 5 MeV, 
and /n|y 30 fii dfi' ~ 0.05. 

Finally, we must estimate the matrix element between the spherical orbitals, 
(iPi,+k\Htp\^2,+k) ■ Using Eq. fll]) for the form of if T p and p(r') = 0(r' — i? )/(| 7rr o 3 ) 
we get 

Htp = ■ n >(^' - *>)■ (38) 

Using Eq. (|3~4"D and (<r • n') 2 = 1 gives 

<^i, + idF T p|^2, + K> = ^ 87rv ^ mro 3 ^i(^o)^(fio)^o 2 « ^1 eV. (39) 



Using \E + - E_\ ~ 50 keV (see, e.g., @,|), & ~ 0.1 (see, e.g., 0), and Eqs. fl35J), 



and (B3) gives (for A « 225) |a| ~ 2f3 3 A- 1 ^r] eV /\E + - £_| ~ 7 x lO" 7 ??. This provides the 



following estimate for the collective EOM: 

|C C oii| ~ 0.05e[3 3 2 ZA 2/3 r 3 r]eV/\E + - E_ \ ~ 4 x lO'V^m) 3 - (40) 

We see that the collective EOM is two orders of magnitude larger than the EOM due to 
unpaired protons. 

We can do a more accurate calculation of the EOM using Refs. |HJ, where the mixing 
coefficients, a for various nuclei were calculated. We use the values from || that were 
calculated using the Woods-Saxon potential. We took the values of the fa parameters from 
fl8| and the nuclei's angular momenta were taken from |I5|]I7|| . The results are shown in 
table [TT] for various nuclei. 

Note the large value of 229 Pa's collective octupole moment. This is due to its large value 
of a, which is caused by the small energy splitting between the members of its parity doublet 
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PH. The possible existence of a static octupole deformation in 229 Pa was stated in J20 . 
However, more recent papers cast doubt on the existence of such a deformation (see, e.g., 
|2T| , |22[] ) . Therefore, the result given for 229 Pa must be understood as being conditional on 



it having a static octupole deformation. 

V. THE ATOMIC ELECTRIC DIPOLE MOMENT INDUCED BY A NUCLEAR 

ELECTRIC OCTUPOLE MOMENT 

In this section we consider the electric dipole moment of an atom that is induced by a 
nuclear electric octupole moment. The electric potential, <f)(r) outside an arbitrary charge 
distribution can be expanded in terms of spherical harmonics as (see, e.g., 



r 



E E ^— • ( 41 ) 

(=0 m=—l 



Here we neglect that part of the potential that comes from the screening of the nucleus's 
Coulomb field by the atomic electrons as it is not important in the octupole potential (see 
the appendix). The qi m are spherical electric multipole moments and they can be written 
in terms of the charge distribution p c (r') as follows [p3|: 



Qim = jY: m {9',<l>y l p c {v')dh'. (42) 
The octupole term of the above electric potential is 

0( 3 )(r) = ^g 3 o^3o(M), (43) 

where we have only taken the m = term, as we will be calculating the matrix element of the 
perturbation between states with the same angular momentum projections. We can write 
this in terms of the octupole moment by using the relation = / Y^O, 0)r 3 p c (r) d 3 r = 
' O zzz . The potential energy, U oc t of an electron (of charge — e) in 0( 3 )(r) will then be 



•1 



U oct = -S^eO±Y ao {0,<l>) (44) 



(for the nuclear I z = I state). 

12 



Using perturbation theory, the electric dipole moment of, for example, an atom with one 
electron above closed subshells induced by U oct can be written as 

7 ,7, ,7, „ (ki\r z \k2)(k 2 \U oc t\h) 

d z = -e(ij\r z \ij) = -2eJ2 ' ' — , (45) 

\k 2 ) ^1 ~ h k 2 

where ip denotes the perturbed atomic wave function, \k±) = |ni, ji, h, m) is the unperturbed 
single-electron ground state, and {I&2)} is the set of states with which \ki) is mixed by the 
perturbation. 

According to the triangle rule for the addition of angular momenta, (ki\r z \k 2 ) can only 
have a nonzero value if \ji — j 2 \ < 1 < ji + j 2 - Similarly, for (k 2 \U oct \ki) to be nonzero, we 
must have \j\ — j 2 \ < 3 < ji + j 2 . This implies that the following conditions need to be 
satisfied for the dipole moment to be nonzero: 

\ji - 32\ < I and jx+j 2 > 3. (46) 

The lowest pair of values that satisfies this condition is j± = 3/2 and j 2 = 3/2. Therefore, 
s states cannot contribute to the dipole moment induced by the nuclear EOM. Also, one 
of l\ and l 2 must be even and the other odd, since the electric dipole moment is a parity 
nonconserving effect. 

We will carry out a relativistic calculation of the matrix element of U oct between the 
single-electron states \ni,jx, li, m) and 1^2,^2; h, m )- (Note that although we wrote Eq. ([45]) 
for an atom with one electron above closed subshells, this single-electron state matrix element 
can be used for all atoms to compare various sources of atomic EDMs.) The relativistic wave 
function of an electron is (see, e.g., |T5| ) 



1 fnfl(r)n jlm (6,<j>) ^ 



njlm 



(47) 



\ 9nji{r)i{-(T ■ n)Qji m (0, 0) j 

where <r is twice the spin operator of this electron and n = r/r. Evaluating the matrix 
element (using (cr • n) 2 = 1) gives 

(nxjilxm\U oct \n 2 j 2 l 2 rn) = -h^eO^hn^Y^j^rr^T, (48) 
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where T is a radial integral: 

r°° 1 2 
^ = Jo ^•4^ ni - ?1 ^ 1 ^ ?n2 -? 2 ' 2 9mjih9najah')'^ ^t. (49) 

Because of the factor of 1/r 4 in the above integral, most of the contribution to T comes 
from small values of r. This allows us to use the following expressions for f n ji and g n ji, for 

fnjl{r) = —[(7 + «)^2 7 (x) - |j 27 -l(x)], 

flWi(r) = — ZaJ 27 (:r), (50) 

where 



!8Zr 



V a 

7 = y/U + \? ~ Z 2 a\ 
«=(-l) i+1/2 -^ + |), 

c -' = rIwJ ' (51) 

where the J's are Bessel functions, a is the Bohr radius, and v is the effective principal 
quantum number {E n \ = —13.6 eV/V 2 ). To avoid confusion, note that I here is the orbital 
angular momentum of the electron, rather than the nucleus, as used in Sec. PH Also, the 
k used here is distinct from the k defined in Eq. flT5p. Carrying out the integration in Eq. 
(f|9|), we obtain 

192(-l)^-* +1 Z 2 r(-3 + 7 i + 7 2 ) , , 

a^vTvT r(4 + 7i + 72)r(4 + 7l - 72 )r(4-7i+7 2 ) 1 ' 

x {(3 + 7i - 72) (3 - 7i + 72) (3 + 71 + 72) (2 + 7i + 72) 

-5(71 + «i)(3 - 71 + 72) (3 + 71 + 72) - 5(72 + « 2 )(3 + 7i - 72) (3 + 71 + 72) 
+30[(7i + /€i)(72 + ^ 2 ) + (^a) 2 ]}, 

where T is the gamma function. 

To illustrate the numerical values involved, we will evaluate a value of the matrix element 
for an electron's interaction with the single particle octupole moment for 209 Bi. This atom 
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has one 6p 3 / 2 electron above closed subshells ((6pi/ 2 ) 2 (6si/ 2 ) 2 . . .)• We will consider mixing 
between d 5 / 2 and p 3 / 2 states, each with an angular momentum projection of 3/2. Using table 
| and Eqs. and © gives 

(n lP3/2 ,m = 3/2\U oct \n 2 d 5/2 ,m = 3/2) « 1.9 x lO- 13 ^)" 3 / 2 ^ -1 . (53) 

We now compare this with the corresponding matrix element for the interaction with a 
nuclear magnetic quadrupole moment (MQM), which is another possible source of an atomic 
EDM. The MQM induced by the T- and P-odd interaction ([!]) was calculated in Ref. ||, as 
well as the matrix element of the interaction between an electron and the MQM field. For 
209 Bi the MQM is « 4.8 x lCr 8 ?7e(fm) 2 . We get the following estimate: 

(n lP3/2 ,m = 3/2\U qnad \n 2 d, /2 ,m = 3/2) « 1.3 x lO^Vi^r^fym -1 . (54) 

By comparing the values of these matrix elements it can be seen that the atomic EDM 
induced by a nuclear EOM will be less than about 1% of that induced by a nuclear magnetic 
quadrupole moment. This same result held for the other atoms and mixing states that we 
considered. Another mechanism which could produce an atomic EDM is the nuclear Schiff 
moment, which for heavy atoms (which we are considering) gives a contribution comparable 
to that of the MQM |J. Because of this we can conclude that the atomic EDM induced by 
a single particle nuclear EOM is negligible in comparison with other possible mechanisms. 
(See the appendix for a discussion of the relative contributions of the octupole, magnetic 
quadrupole, and Schiff moments.) 

Now we will give example values of the matrix elements for those nuclei having a static 
octupole deformation. As well as possibly having collective EOMs (as shown in Sec. |TV| ), 
these nuclei can also have collective MQMs that are an order of magnitude larger than the 
single particle MQMs discussed above ||24|| . We will use 225 Ac as an example and once again 
we will consider mixing between d$/ 2 and p 3 / 2 states, each having an angular momentum 
projection of 3/2. This nucleus has an MQM of ~ 2 x 10 _7 r]e(fm) 2 |2"4] . We obtain 

\( ni d 5/2 ,m = 3/2\U oct \n 2 p 3/2 ,m = 3/2) | « 1.4 x lO" 11 ^)"^^- 1 (55) 
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and 

|(n 1 tZ 5 /2 > m= 3/2|f/ quad |n 2 p 3 / 2 ,m = 3/2) | ~ 6 x lO^^^-^cm- 1 . (56) 

These results show that the contribution of the collective nuclear EOM to the atomic EDM 
is smaller than that of the collective MQM. This also applies to the other isotopes shown in 
table HI, except for 229 Pa. ( 223 Rn does have a fairly large collective octupole moment, but 
since it has a closed electron shell the EOM does not contribute to the atomic EDM.) For 
229 Pa the contribution of the EOM may be comparable to that of the MQM, but as stated 
above, it is not certain that this nucleus has a static octupole deformation. 

VI. POSSIBLY ENHANCED MAGNETIC QUADRUPOLE MOMENTS IN 
NUCLEI WITH AN OCTUPOLE DEFORMATION 

Finally, we discuss a mechanism by which single particle magnetic quadrupole moments 
could be enhanced in even-odd nuclei with an octupole deformation. We will first consider 
the MQM of such a nucleus in its intrinsic (body-fixed) frame and then transform the MQM 
into its laboratory frame. 

As in Sec. [IV] the wave function of the external nucleon in the intrinsic frame is ^^(r'), 
defined by Eqs. ( |32|) , (|33|) , and fl34|). The MQM in the intrinsic frame is then 

M intr = (i} +K \M z > zl \^ +K ) (57) 
= 2j{ik t+K \M z . gl \ifo t+K ), (58) 

where M is the operator for the MQM (see Ref. 0). Note that the MQM is defined for the 
maximum value of the projection of the angular momentum onto the z'-axis, so K = j = I 
(for the ground state of the rotational band I = j). Using a result from Ref. || we have 

M ntr = 7^r— 7 1 -<> - / R i R * r r2dr > ( 59 ) 
1 + 1 m J 

where \l is the magnetic moment of the external nucleon in nuclear magnetons and q = (1) 
for a neutron (proton). It should be noted that this (intrinsic frame) single particle MQM 
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differs from the spherical nucleus single particle EOMs and MQMs (considered in Sec. p 1 1 
and Ref. |J, respectively) in that the former is generated due to the interaction V3, coming 
from the nucleus's octupole deformed shape [see Eqs. ([33]) and ([36])], while the latter is due 
to the interaction Htp ([!]). The Htp interaction will come into the current situation when 
we transform into the lab. frame. 

Since we are only working out a rough estimate of the MQM we take the integral in the 
above equation to be 

c 1 

J R 1 R 2 rr 2 dr ~ -r A 1/3 , (60) 

where r$ = 1.1 fm. For both types of nucleons we have \fi — q\ ~ 1.8 and so we can write 

|M intr | ~ .2A 1 /3^_l| 7 |e(fm) 2 . (61) 

Now we turn to the MQM in the laboratory frame. The wave function of the nucleus 
in the lab. frame is (as in Sec. |TV| ) if) = ip + + atp~ , where ip^ is defined by Eqs. ( p5|) and 
(fZ3|) and a is given by Eq. Q2"7|) — this is where the interaction Htp comes into the present 
situation. The MQM in the lab. frame is then 

M lab = (iP\M zz \iP) 

= a((IMK\M zz \IMK) - (IM - K\M ZZ \IM - K)) 

= 2a(IMK\M zz \IMK). (62) 

(Once again, the MQM is defined for M, the projection of the angular momentum onto 
the z-axis, equal to I.) Note that (IM - K\M ZZ \IM - K) = -(IMK\M ZZ \IMK), as 
ip±K = ipi,±K ±^ip2,±K and so (i)- K \M z ' z >\ip-. K ) = -(i) +K \M z > z '\i/} +K ). 
Now consider (IMK\M ZZ \IMK) . We have [from Eq. @§] 

(IMK\M ZZ \IMK) = J D^ K (<p,e,0)r K (r')M zz (e, ^0)^)01^^,6,0) d 3 r' dQ. 

(63) 

17 



Transforming M zz from the lab. frame to the intrinsic (x',y',z') frame gives M zz (6,<ft,0) 
Dl*((p, 9, 0)M z , z r(9', </>', 0) (see, e.g., [fL4| ). Substituting this into Eq. (H) and using Eq. (| 



gives (using a formula for the integral of D-functions from [[14]]) 



(IMK\M ZZ \IMK) = 2 I±±M- mtT I D^ K ( ( f>,e,o)D 2 M,9,o)D I MK ( ( p,e,o)dn 



4tt 

M intr (J, M- 2, 0| J, M) (I, K- 2, 0| J, X) 

\2 



M- mtr {I,I;2,0\I,I) 
1(21 - 1) 



(J + l)(2/ + 3) 
(Recall that M = K = I.) Using Eq. (H) then gives 



M intr . (64) 



Mlab = 2a (/ + ^ + 3) Mto (65) 



and so [using Eq. (|6l"l) 1 



|M lab | ~ O.^jalM (j+ip^ + g) 6 ^) 2 - (66) 



Now we have |a| ~ 7 x 10 7 r] and [7) ~ 0.1 (see Sec. [IV|) . For A « 225 we have 



|Mi ab | ~6x 10~V(ftri) 2 - (67) 

This is smaller than the collective MQM due to the spin hedgehog mechanism — ~ 2 x 
10 _7 ?7e(fm) 2 |24| and so it is not the dominant mechanism. It is, in fact, of the same order 
of magnitude as the single particle MQM [e.g., for 209 Bi, approximately 4.8 x 10 _8 ^e(fm) 2 
(see Sec. 0)]. This means that there is no enhancement. 

Enhancement was a possibility here due to the relatively large value of a that comes from 
the smallness of the energy splitting between members of the parity doublet (E + — EJ). 
However, the inclusion of the factors of [7! ~ 0.1 and 1(21 - I) /[(I + 1)(2J + 3)] ~ 0.3 
(this enters on transforming to the lab. frame) ensures that the possible enhancement is not 
realised. 

However, 229 Pa may be an exception to this, as it has a large value of a, but it must be 
remembered that it is not certain that this nucleus has a static octupole deformation (see 



Sec. (TVI) . If is does have such a deformation then its MQM due to the present mechanism 
would be |M lab | ~ 3 x 10~ 7 ^e(fm) 2 [using Eq. (0), the value of a given in Ref. 0, and 
|7| ~ 0.1], which is of the same order of magnitude as the collective MQM due to the spin 
hedgehog mechanism. 

Note that the enhancement of MQMs in deformed nuclei with opposite parity levels close 
to each other has also been considered in Refs. {§J|. 
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APPENDIX A: THE RELATIVE CONTRIBUTIONS OF VARIOUS NUCLEAR 
MOMENTS TO THE ATOMIC ELECTRIC DIPOLE MOMENT 

An atomic electric dipole moment can be generated by T-,P-odd nuclear moments. How- 
ever, according to the Purcell-Ramsey-Schiff theorem |25],|26|, a nuclear EDM cannot generate 
an atomic EDM. This happens because the electrostatic potential of a nucleus is screened 
by its atomic electrons. Therefore we must look to higher order moments. 

The (screened) electrostatic potential of the nucleus can be written as (see, e.g., H@,H) 

where p(r) is the nuclear charge density (/ p(r) d 3 r = Z), d = J erp(r) d 3 r is the nuclear 
EDM, and Vj = di = The first non-zero T-,P-odd term in this potential is |||8],|9| 



\(3) 



lr 1 1 r 1 

-- J ep(Y)r i r j r k d 3 rdidjd k — + — J p(r)rirj d 3 r (d • V)^ — 

1 r 1 1 

-- J ep(r)\r i r j r k - -r 2 (riS jk + Tjb ik + r k Sij)) d 3 rdidjd k — 



R 
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If 1 

+ ^Jp(r)r%d 3 r(d-V)d l d~ 
~Oi jk didjd k - + ~^Qa ( d • V)didj^ 

-— J ep(r)r 2 ri d 3 r diS jk djd k - + — J p(r)r 2 d 3 r (d • V)^-^^— 
^octupole + ^Schiff' ( A2 ) 



where 



d>^ \ i — — OjikdjOnOb— H ^QadhOjOnOb — 

^octupole q 3 3 ji e2Z R, 

~ -\o ijk didjd k —, (A3) 

4 3 C hiff = - s • vy2 ^ = 4?rS ■ v w ( A4 ) 

(using SjkdjdkR' 1 = V 2 i? _1 = — 4tt5(R)). Qij is the T,P-even electric quadrupole moment: 
Qij — I e p( r )( r i r j ~ \ r2 &ij) ^ 3r - The vector S is the nuclear Schiff moment and the rank 3 
tensor Oy/. is the nuclear electric octupole moment. These are given by 

Oijk = J ep(r)[rirjr k - \r 2 (ri5 jk + rj5 ik + r fc <%)] d 3 r, (A5) 

S = — ( I ep(r)r 2 rd 3 r - - A— f p{v)r 2 d 3 r) . (A6) 
10 \J 3 Z J J 

The second term in the octupole potential (|A3|) comes from the screening of the nucleus's 
Coulomb field. Because only the non-spherically symmetric part of the density will give a 
non-zero value of Qij, only the external proton contributes to this screening term (in the 
case of a deformed nucleus all the protons in the external shell contribute). This, together 
with the presence of the factor 1/Z, means that it will be small and so we neglect it. The 
screening term in the Schiff moment [the second term of Eq. (|A6|) 1 is not negligible as the 
whole density p(r) contributes to it and so all of the protons are involved. 

^Schiff anc ^ ^octupole ^°^ n a PP ear a ^ the same order in the expansion of Eq. (|Alj), i.e., 
both potentials contain third derivatives of 1/R and their corresponding nuclear moments 
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are both integrals with integrands containing r 3 terms. Therefore, we should expect the 
contributions of the Schiff and octupole moments to the atomic EDM to be roughly the 
same. However, since the octupole moment is a rank 3 tensor, as opposed to the (vector) 
Schiff moment, an electron interacting with an octupole moment must have a higher angular 
momentum (see Sec. [V]) than it would need for a Schiff moment. This results in the octupole 
moment's contribution to the atomic EDM being smaller than the Schiff moment's, as a 
higher angular momentum wave function penetrates the region close to the nucleus less due 
to the greater centrifugal barrier. This is confirmed by the calculations in this paper (see 
Sec. 0). 

An atomic EDM can also be caused by a nuclear magnetic quadrupole moment (MQM). 
The vector potential due to such a MQM is 

Ai(R) = -^iiA^i (A7) 
o R 

where M kn is the (rank 2) magnetic quadrupole moment tensor: 

M kn = - J {r k e npq + r n e kpq )j p r q d 3 r (A8) 

(j p is the electromagnetic current.) Note that this is a lower order term than the Schiff and 
octupole moments — the potential contains second derivatives oil/ R and the integrand in 
Eq. (|A8|) contains r 2 terms. However, the contribution of the MQM to the atomic EDM 
will actually be of roughly the same order of magnitude as that of the Schiff and octupole 
moments. To see this, consider the (nonrelativistic) interaction between the atomic electron's 
magnetic dipole moment and the magnetic field from the MQM: fi e B = /i e (V x A), where 
He is the Bohr magneton. Now from Eq. (TA7|) we have V x A ~ V x (M/R 3 ) ~ M/R 4 , 
where M is the magnetic quadrupole moment. From Ref. || we have M ~ Hn£, where /xjv 
is the nuclear magneton and £ is defined by Eq. (|3]). Therefore we have (in h — c — 1 units) 

[i e [i N eel 
^ B ~^~-- W t (A9) 

Now we compare this with the electron's interaction with the octupole: 
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^octupole ~ & ~ -rT-^' ( A1 °) 

using the result for the octupole moment in Eq. (pi]) (i?jv = 1.L4 1 ' 3 fm is the radius of 
the nucleus). Therefore, the ratio of the contribution of the octupole moment to the atomic 
EDM to the contribution of the MQM is ~ R^ 2 m e m p ~ 10~ 2 A 2 / 3 . A more accurate estimate 
that takes into account the electron angular momenta dependence of the matrix elements 
gives a ratio that is closer to 10~ 4 A 2 ^ 3 . In Sec. the ratio was ~ 0.01, which is consistent 
with this result. 
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TABLES 



Nucleus Proton state O sing [^(fm) 3 ] 

209 Bi hg /2 1.2 x 1(T 7 

133 Cs g 7/2 7 x 1(T 8 

127 I d 5/2 -1.3 x 1(T 7 



TABLE I. Approximate values of the unpaired proton octupole moment (O s ing) for various 
nuclei with I > 3/2, in terms of the parameter 77. 



Nucleus / I e> coU | [^(fm) 3 ] 

223 Ra § 2 x 1(T 6 

223 Rn 1 7 x 1(T 5 

225 Ac § 7 x 10- 6 

223 Fr 3 4 x 1Q -6 

229 Pa I 3 x 10- 4 



TABLE II. Approximate values of the collective octupole moment (O co n) for various nuclei 
that have a static octupole deformation and / > 3/2, in terms of the parameter 77. 
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